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Renormalization and Universality of Van der Waals forces* 
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Abstract. Renormalization ideas can profitably be exploited in conjunction with the superposition principle of 
boundary conditions in the description of model independent and universal scaling features of the singular and 
long range Van der Waals force between neutral atoms. The dominance of the leading power law is highlighted 
both in the scattering as well as in the bound state problem. The role of off-shell two-body unitarity and causality 
within the Effective Field Theory framework on the light of universality and scaling at low energies is analyzed. 
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1 Introduction 

Van der Waals (VdW) forces were first conjectured from 
the experimental observation that in an adiabatic expan- 
sion a gas of neutral particles cools down (Joule-Thomson 
effect). Since the inter-particle distance at room temper- 
ature is ~ 30A this suggests that VdW forces are long 
range and attractive. Q Their genuine quantum mechani- 
cal origin and form ~ l/r^ was uncovered by London ||2l 
as long range dipole fluctuations between charge-neutral 
atomic and molecular systems. They dominate at distances 
above 5 - lOA and hold atomic dimers together The rel- 
flf/vwf/c Casimir-Polder forces ~ 1/r^ include retardation, 
are a consequence of vacuum fluctuations 1 3 1 and operate 
at very long distances ~ 1000 - 2000A, a relevant scale in 
colloids. The general field theoretical treatment due to two 
photon exchange [41 yielded the so far missing magnetic 
contribution (For a review see e.g. |5|). 

Van der Waals forces, besides being long range, di- 
verge if directly extrapolated to short distance scales but a 
sensible interpretation becomes possible [6 7|. Because of 
the interest on ultra-cold atoms in recent years |81 funda- 
mental work for neutral atoms was initiated in Refs. 0911011 II 
(see also lfT2l ') from the point of view of quantum defect 
theory, where a spectacular reduction of parameters takes 
place. This is supported by more conventional potential 
calculations and a pattern of (VdW) universality and scal- 
ing sets in 1 1 3 1 with no explicit reference to short distance 
scales or cut-offs. Actually, Effective Field Theories (EFT) 
explicitly exploit the characteristic low energy parameter 
reduction from the start and yield very general universality 
patterns which do not resolve the nature of the forces and 
therefore enjoy a wide applicability 0141151161 . They are 
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based on pure contact (zero range) interactions and discard 
the long distance tail of VdW forces. In the present contri- 
bution we analyze the quantum mechanical problem from 
the point of view of renormalization, and address to what 
extent do these contact interactions faithfully describe the 
underlying Van der Waals force. 



2 From Binding to Van der Waals forces 

To provide a proper perspective it is interesting to recall 
the distance scales where we expect the dispersion forces 
to dominate. For simplicity let us consider the simplest 
H2 molecule, which Hamiltonian in the CM frame and 
in the Born-Oppenheimer approximation, valid for heavy 
protons, rup » m^, reads 

H^Hi+H2 + Vn (1) 
where the single atom hydrogen-like Hamiltonians are 



Hi,2 - --. 



1 



- 

1,2 



2m, '-^ |ri,2 + r/2| 
with ri and the electron coordinates and 

1 1 



V17 - a 
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(2) 



(3) 



Defining ri -t = ri+r/2 andr2,± = r2±r/2, the solutions to 
Eq. (|2]i are i/'„(ri__) and 4'm(X2,+) where £„ - —m^a^ I2n^ - 
-13.6eV/n^ [17] so that for Vn - Q the total mirror sym- 
metric molecular normalized wave function reads 



iAn(ri -)i/r,„(r2,+) ± i/',„(ri,+)tA„(r2 -) 
V2(l ±5„„,) 



(4) 



^ We work in natural units with % = c = \ and a = 
liAneohc) = a = 1/137.04 the fine structure constant and 
he = 1973.2AeV. The Bohr radius is a,, = n^/(m,a) = 0.51A. 
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Fig. 1. Bom-Oppenheimer interaction energies for thie H2 
molecule (in units of the ionization energy of the H atom, = 
-nigC^a-/! = 13.54eV) as a function of distance (in units of the 
Bohr radius oo = ti^/m^c^ = 0.5 1 A). We compare first order, zliii , 
and second order, AE2, perturbation theory with the pure VdW 
approximation zliivdw = -6mgC^a^(ao/r)^. All calculations use 
the closure approximation. 



with S„m(r) the corresponding overlap integral, fulfilling 
Sn,miO) = <5n,m and S „^,„ - 0{e^^''^"°). This generates a 
coupled channel matrix Hamiltonian what eigenvalues pro- 
vide the H-H adiabatic energy, £'hh('") — » 2£'h for r — > 00. 
The nowadays standard variational approach pioneered by 
Heitler and London |18 | and culminating with the bench- 
mark determination of the ground state dissociation en- 
ergy 1 19 1 does not accurately work at very long distances, 
and thus perturbation theory might be preferable. Taking 
//() - Hi + H2 as the unperturbed Hamiltonian and Vn 
as the perturbation, one can determine the potential energy 
shift of the system at a fixed proton-proton separation r in 
perturbation theory, which to second order reads, 

VH-H(r) = Eii-Hir) - 2Eh = zl£i(r) + mir) + ... (5) 

In the case of the H2 molecule the calculation was under- 
taken in 1930 by London and Eisenschitz |20| in the clo- 
sure approximation (CA) 0. We have reproduced the an- 
alytical calculation ETll and the results are presented in 
Fig. [U Already in zfisiCr) the finite atomic size yields ef- 
fects which are ~ g^^r/oo 

An interesting feature is that due to the finite atomic 
size, ~ fl(), second order in perturbation theory is finite at 
zero separation, |zl£'2(0)| < 00. Actually, for r — > we 
expect the exact behaviour Vu-Hif) - a/r — > Eue - 2£'h- 
In the CA we get E'j^^ - -llOeV to be compared with 
= -79eV. At large distances only the second order 



^ This somewhat crude approximation corresponds to replace 

y JV^^(y^)os>2K> 
~^ E„ - Eq El - Eq 

Note that the sum includes also continuum e - p states. 



direct term in Eq. (|5]l contributes yielding, 

V(r) = -2-f +0(e-2^/"») (6) 

n=6 

where in the CA C„ = -c„meC^a^a".^ with ce = 6, eg = 135, 
cio = 7875/2, about 5% accurate compared to exact val- 
ues (see table [Til. For n » 6 one has C„ ~ (ao)"nl which 
means that the VdW series represents a diverging asymp- 
totic expansion. As it usually happens in such a case, for 
a fixed order of the truncation, this poses a lower limit on 
the distance below which it makes no sense improve the 
calculation. Usually, only the terms with n = 6,8, 10 are 
retained. Calculations show that, so far, corrections are 
always negative, i.e. C„ > as could be infeiTed from the 
Lieb-ThiiTing universal bound |23| which establishes that 
generally V(r) < -c/ for any pair of atoms with under- 
lying Coulomb forces such as in Eq. The positivity of 
the C„'s will prove essential in what follows. As we see in 
Fig. [T]in the H2 case the Van der Waals force dominates 
for distances of about 5A. Actually this corresponds to an 
energy which is comparable to the environmental transla- 
tional thermal energy kgT ~ 1/40 eV [24 1. Clearly in the 
ultra-cold region the VdW force dominates. Therefore, a 
VdW theory should work in a very broad range of dis- 
tances and wavelengths. 

The VdW potentials are valid assuming large distances 
r » ao and low energies p^/(2jj) A - Ei - E^ so that 
intermediate excited states do not contribute dynamically. 
Thinking of the molecular wave function like e.g. Eq (|4|l, 
at large distances we may neglect the exchange term with 
accuracy (9(e ''^''°), and the coupled channel Hamiltonian 
reads 

This problem may quite generally be decomposed within 
the total multichannel Hilbert space into the P-space (elas- 
tic) and g-space (excited) states [25 1, where P - |0, 0)(0, 0| 
and Q - I - P are the corresponding projection operators. 
This yields the box-matrix structure 

HppWp + VpqWq^E^p 
VqpWp+HqqWq^EWq (8) 

where in the Born-Oppenheimer approximation all these 
potentials are local functions of the inter-nuclear separa- 
tion r. Eliminating the unobserved channels we get the ef- 
fective optical potential 

Vpp{E) = Vpp + Vpq{E - HqqY'Vqp (9) 

which develops an imaginary part if the first inelastic thresh- 
old becomes open. The important point is that if the com- 
plete underlying electronic Hamiltonian, Eq. ([T]i is energy 
independent, then for E - p^/(2jj) + 2Eu < Ew + Eu 

V'pp(E) = -Vpq{E - HqqY^Vqp < (10) 

See e.g. the impressive calculation in Hydrogen up to 
C32 1221 for which we can fit a behaviour C„ ~ (l/2)"«! in atomic 
units. In the CA we have = 1.93 x 10"^ to be compared with 
Cg""^' = 2.51 X 10^* Ga. 



19 International lUPAP Conference on Few-Body Problems in Physics 



The underlying local and energy independent dynamics 
has consequences in the low energy EFT representations 
of the VdW interactions (see Section|6]l. 



3 Universal Scaling Theorems 



We review here some results found in a previous work 1261271 
(for a short review see e.g. Il28l ).within a nuclear physics 
and multichannel context which will prove useful in the 
analysis of VdW forces . Our starting point is the finite en- 
ergy scattering state Schrodinger's equation for the relative 
wave function between two particles of masses mi and ma 
which interact through a central potential, 



1(1 + 1) 
• u + U(r)u H T — u 



(11) 



where U(r) - 2fiV(r), fi - m\m2l{mi + 11x2) is the reduced 
mass, k - plfi - InjA the wavenumber and u{r) the re- 
duced wave function. We will neglect finite size and ex- 
change effects and take V{r) given just by Eq. (|6]l where 
C„ are the standard VdW coefficients which are computed 
ab initio from electronic orbital atomic structure calcula- 
tions [29 1 . The potential in Eq. ^ can conveniently be 
rewritten as 



Uir) = —I 



1 + gx — + 02— + ■ 



(12) 



where /?6 - (2yuC6)^ is the VdW length scale and g\, §2, 
etc. represent the contributions from Cg, Cio etc. at r = 7?6 
respectively. In table[T]we compile numbers for a bunch of 
interesting cases. Typically, R(, ~ \Q - 200A but also g\ ~ 
10"^ and 02 ~ 10"'*. This raises immediately the question 
under what conditions the expansion (fT2] i can be truncated 
in a meaningful way, i.e. when the neglected terms can 
indeed be considered negligible in scattering and bound 
state properties. This question is intriguing since at short 
distances the more singular terms are manifestly more di- 
vergent. Actually, the range where Cg yields an important 
correction but Cio is still small is g^^^* <k r/R^ <K g^^^^ 
which in view of table [T] g2 ~ g\, becomes extremely nar- 
row or inexistent. 

To determine the solution of Eq. ( fTTT i it is necessary to 
give sensible boundary conditions at the origin and infin- 
ity. For the usual regular potentials there are a regular and 
irregular solution at the origin, and the regularity condition 
m(0) = fixes uniquely the solution. However, since the 
potential is singular and attractive there are two linearly 
independent solutions, so regularity at the origin does not 
select a unique solution. Indeed, at short distances the De 
Broglie wavelength is slowly varying, d[U{r)]^^ /dr <s: 1 
and hence a WKB approximation holds {6'T\, yielding for 
r ^ 



Uk(r) ^ C — 



n/4 



sm 



2 IR„\ 
n-2\ r I 



(13) 



where R,, = (l/uCny^^" corresponds to the highest VdW 
scale considered in Eq. ^ (see also Eq. (fT2l i. The phase 



ifik is arbitrary and could, in principle, be energy dependent 
(see below. 

To fix ideas we will restrict to / = 0, i-waves. For a 
positive energy scattering state it is convenient to use the 
normalization at very long distances given by 

smikr + 60) n \ . I .X sin(fcr) 

Ukir) — » = cos(kr) + k cot do , (14) 

sm<5o k 

where Soik) is the scattering phase shift for the / - an- 
gular momentum state. For the potential which at long dis- 
tances behave as Eq. (fT2l i one has the effective range ex- 
pansion (ERE) 1361 

kcotSoik) = -— -H -rok^ + V2k'^log{k^) + ... (15) 
ao 2 

where q-q is the scattering length, and ro is the effective 
range, q-q and ro can be calculated from the asymptotic be- 
haviour of the zero energy solution 



uoir) 1 - r/ao , 
and using the definition 



(16) 



dr[il-r/aof-uo(rf]. (17) 



Next, we use the superposition principle of boundary con- 
ditions and write 



Ukir) = M/t,c(r) -I- k cot 60 Ukjf) , 



(18) 



with Uk^cif) — > cos(A;r) and Uk.s{r) — > sin(kr)/k for r — > 00. 
At zero energy we have 



uo(r) = Mo_,.(r) - U(iJr)/ao , 



(19) 



with Mo,e('") 1 and mo,j('") — > sinr for r ^ 00. The short 
distance phase ipo can be fixed in practice introducing a 
short distance cut-off, re- The way to proceed is as follows. 
Given a scattering length ao as input one integrates in from 
large down to small distances, say r - R^ whence 

determining tpo. To determine tpk we do the same but for 
finite energy states. A relation between both short distance 
phases can be found as follows [26). If we build (m^mq - 
MqMo)' and integrate from to infinity and use Eq. ( fTsT i 
and Eq. (T9[ respectively, one gets 



^ siniipk - ipo) 



^k^ \ dr 
Jo 

y[ukAr) + k cot Soik) Uk/r)\ , (20) 



M(),c('") uo^sir) 



which becomes an orthogonality relation if and only if ipk - 
ifio- This implies that the corresponding Hamiltonian, al- 
though unbounded from below becomes self-adjoint on the 
domain of square integrable functions which have the short 
distance behaviour, Eq. ( fTST l. with the same short distance 
phase (a common domain of definition). In other words, the 
short distance common phase, ^, labels the particular self- 
adjoint extension of the Hamiltonian depending on the Van 
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Table 1. Reduced di-atomic masses, Van der Waals dispersion coefficients Ce, Cg and Cio as well as the leading Van der Waals length 
scale R(, = {2jiC(,)^ , and the dimensionless coefficients gi and defined by the dimensionless reduced potential 2fiV(r)R^ = R^U(r) = 
-x"^ 1^1 + gix^^ + 92 * + ■ • ■ ] with X = r/Rf, the distance in Van der Waals units, of the di-atomic systems used in the present pa- 
per. Atomic units are used throughout. The atomic masses have been taken from the National Institute of Standards and Technology 
[http://physics.nist .gov/PhysRefData/Compositions/' 



Atoms 






C8(10'a.;/.) 


Clo(10^(3.M.) 




gi(10"^) 


92(10"'*) 


H-H 


918.576 


6.499 130] 


0.001244 |30| 


0.0003286 f30l 


10.4532 


17.51760 


423.45441 


He-He 


3648.150 


1.461 |30| 


0.000141 t30J 


0.00001837 |3_0J 


10.1610 


9.35937 


117.94642 


Li-Li 


6394.697 


1389. 1311 


0.834 (32) 


0.735 |32j 


64.9214 


1.42458 


2.97874 


Na-Na 


20953.894 


1556. |31| 


1.160 (32] 


1.130 1321 


89.8620 


0.92320 


1.11369 


J\.-JV 


JJJLJ.Z.H- 1 


Joy 1 . 




^ ^70 




AzL7S0 


ZLQ78A 


Rb-Rb 


77392.363 


4691. |31| 


5.770 1132] 


7.960 |32l 


164.1528 


0.45647 


0.23370 


Cs-Cs 


121135.907 


6851. LIU 


10.200 (32) 


15.900 |32| 


201.8432 


0.36544 


0.13983 


Fr-Fr 


203270.053 


5256. (21] 


6.648 133] 


10.699 |33| 


215.0006 


0.27362 


0.09526 


Li-Na 


9798.954 


1467. (31] 


0.988 |32] 


0.916 |32| 


73.2251 


1.25605 


2.17183 


Li-K 


10837.871 


2322. (31] 


1.950 (32] 


2.100 | 32 | 


84.2285 


1.18374 


1.79689 


Li-Rb 


11813.296 


2545. p3ll 


2.340 (32] 


2.6101321 


88.0587 


1.18572 


1.70555 


Li-Cs 


12148.102 


3065. 1311 


3.210 (32] 


3.840 |32| 


92.8950 


1.21364 


1.68241 


Na-K 


26356.596 


2447. 131] 


2.240 132] 


2.530 |32| 


106.5708 


0.80600 


0.80155 


Na-Rb 


32978.812 


2683. Ii31j| 


2.660 132] 


3.130132] 


115.3377 


0.74528 


0.65923 


Na-Cs 


35727.672 


3227. (31] 


3.620 (32] 


4.5501321 


123.2277 


0.73874 


0.61148 


K-Rb 


48685.873 


4274. (31] 


4.930 (32] 


6.600132] 


142.8292 


0.56543 


0.37106 


K-Cs 


54924.387 


5159. (31] 


6.620 (32] 


9.400 [32j 


154.2909 


0.53903 


0.32152 


Rb-Cs 


94444.928 


5663. 1311 


7.690 f32] 


11.300 1321 


180.8480 


0.41520 


0.18654 


Cr-Cr 


47340.881 


733. L35J 


0.750 135] 




91.2731 


1.22821 





der Waals dispersion coefficients. 0. One important prop- 
erty is that (fi is not only energy independent but it must be 
fixed independently on the long distance potential, and thus 
encodes all short distance information not accounted for by 
the truncated potential Eq. ( fTSl i. Moreover, there is a fur- 
ther remarkable consequence of the energy independence 
of (f obtained by expanding the integrand in Eq. ( |20] |. 



A:cot<5o(A:) 



aoC(k) + D(k) ' 



(21) 



whereas the functions yi, S, C and D are even functions of 
k which depend only on the potential and are given by 



C(k) 
D(k) 



y-»00 

I dr UQ^c(r)uk,c(r) , 
Jo 

I dr UQjr)uk^c(r) , 
Jo 

/-»oo 

I dr UQ/r)uk,c(r) , 
Jo 

dr uo^s(r)ukjr) . 

a 



(22) 
(23) 
(24) 
(25) 



Note that the dependence of the phase-shift on the scat- 
tering length is wholly explicit; cot So is a bilinear rational 
mapping of q-q- This is just a manifestation of the under- 
lying Moebius transformation well known from the theory 
of ordinary differential equations discussed in lIZTl . The 



' The energy independence can also be deduced from the 
smallness of the wave function at small distances. 



obvious conditions Jl{Q) = D{0) = and S(0) = C(0) = 1 
are satisfied. The limiting procedure poses no problem in 
principle, and the divergence of the potential at short dis- 
tances has been eliminated by demanding a finite physical 
scattering length. This is equivalent to a renormalization 
condition at zero energy. We will analyze below other al- 
ternative renormalization conditions based on one bound 
state. The effective range is defined by Eq. ( fTTl i. and using 
Eq. (1% we get 



B 



C 



where 



ro^A + — + — , 

ao ai 



A = 2 c/r(l-M^^), 
Jo 



B 



-4 dr(r - uo,cUo.s) , 



'I 

Jo 



C = 2 dr(r^-ulJ, 



(26) 

(27) 
(28) 
(29) 



depend on the potential parameters only. The interesting 
feature of the previous equations is that all explicit depen- 
dence on the scattering length oq is displayed by Eq. (l26l). 
This is a universal form of a low energy theorem, which ap- 
plies to any potential regular or singular at the origin which 
falls off faster than l/r^ at large distances. Since the poten- 
tial is known accurately at long distances we can visualize 
Eq. (|26] | as a long distance (VdW) correlation between ro 
and ao. 
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The above results for the effective range and phase- 
shift, Eq. (|26] | and Eq. $2T[ are completely general. If, in 
addition, the reduced potential depends on a single scale 
R, i.e. U(r) = -F(r/R)/R^, one gets universal scaling re- 
lations 



ro - - R -R^ 



R 



ao 



and 



RkcotSo{kR,ao/R) = 



aoJl(kR)+R!B(kR) 
aoC{kR) + RDikR) ' 



(30) 



(31) 



where A,B and C are purely geometric numbers, and ,S, 
C and D are functions which depend solely on the func- 
tional form of the potential. Thus, if a potential has a sin- 
gle scale the phase-shift can be computed once and forever 
for a given energy. This allows to compare quite different 
physical systems which have identical long range forces 
but different short distance dynamics. A remarkable con- 
sequence of Eq. ( l30t is that if ao » R then ro ~ R whereas 
ao <^ R implies ro » R. 



4 The Pure VdW force 

For the pure l/r^ case one has an analytical solution for the 
zero energy state so that the effective range can be com- 
puted analytically fllOlllll using Eq. ( flTl l to get 



37T 



B 



-, C 



37T 



(32) 



in agreement with the general result, Eq. ( |30] |. In Fig. |2] 
we confront the prediction for the effective range to the 
result of many ab initio calculations. As we see the agree- 
ment is rather impressive, meaning that for many practical 
purposes Eq. (l32l i and Eq. (l30l l summarize the relevant in- 
formation on the, about a hundred, data. In Fig.|3]we show 
the VdW universal functions yi(/t), !B(k), C(k) and Dik) 
uniquely determined by the VdW potential and which need 
that the scattering length be specified separately to obtain 
the phase shift (see Eq. dSTT i. 

One important issue has to do with the cut-off depen- 
dence. From a Callan-Symanzik renormalization group type 
argument one has ll27l 



d6(k, rc) _ _ ^3 / "k(rc, rc) 
drr 



(33) 



when ao is fixed. This in turn means that for U{r) R^Jr^ 

and using Uk{rc,rc)/k ~ (rc/R(,)^^^/kone has that finite cut- 
off corrections scale as 



6(k)-6{k,r,)^0{k'r^./R6) 



(34) 



Thus, the short distance cut-off r^ is a parameter for kR(, » 
1 but becomes innocuous otherwise when r^ <k Re, imply- 
ing that universality is robust. We show for a sample case in 



Pi 





Fig. 2. The VdW universal effective range in units of the VdW 
radius R defined by 2/jV(r) = -R*/r^. Compared to several cal- 
culations |37J (Li-Li,Na-Na), |38| (Na-Na), |39| (Cs-Cs), |40| 
(Cs-Cs), ED (Na-Rb), (42j (Be-Be), |43 | (Cs-Rb), |44 | (Cr-Cr), 
I45l (Fr-Fr), (H-H), |?7J (H-H), t48J (H-H), t49J (H-H). 

Fig. |4]the phase shifts for fixed energies, where the rather 
smooth converging pattern can be clearly confirmed. 

Turning to the negative energy bound states with E = 
—Tp-y^ 1(2^), the behaviour at long distances is 



Uy{r) Ae 



(35) 



whereas at short distances the energy dependent boundary 
condition, Eq. (fTSl l , holds. Orthogonality to the zero en- 
ergy state requires that 



ao = 



druy(r)uojr) 
Jq dr Uy(r)uo,c{r) 



(36) 



This relation explicitly yields the scattering length from a 
given bound state. Inverting the formula we get, the energy 
spectrum, E„ - -tp-yljlji for a fixed value of the scatter- 
ing length and a fixed potential^. This relation determines 



^ Using the WKB method one obtains (see also Ref. |27|) 
3^(1) 



n + c 



1 



{jnRV = 0.20587(r„i?)5 



where c is a constant of order unity. Thus if we happen to have 
a zero energy state, 70 = then c = and the bound states are 
E„ = - 1 16.08n' Ti^ liljiR^) counting downward in the spectrum. 
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Universal function A(i;) Universal function B(k) 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

Universal function C(k) Universal function D(k) 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

kRj kRs 



Fig. 3. The VdW universal functions in units of thie VdW radius 
R defined by 2fiV(r) = -R*/r^. Using these functions one can 
determine the phase shift if the scattering length is known by the 
formula A: cot (5o(A:) = [aoJi(kR) + S(kR)]/[aoC(kR) + D{kR)]. 



























kR,= 1 — 






kR6= 5 






kR5= 10 






kRj = 20 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Fig. 4. Convergence of the phase shift for a fixed momentum k as 
a function of the short distance cut-off, in VdW units. 

the energy spectrum E„ from the scattering length ao and 
the potential. Conversely, we can get the scattering length 
ao from a given bound state wave number y. In Fig.|5]we 
display the lowest bound states as functions of oq for the 
pure VdW potential in scaled units. Such an universal plot 
allows also to deduce the scattering length from the knowl- 
edge of the weakly bound states in a complete model inde- 
pendent way. 

We may ask how should the neglected higher order 
1/r^ ,1/?"'", etc. terms in the potential affect the lowest or- 
der 1/r^ calculation. To this end we show in Fig. |6] the 
effect of adding the term 1 jr^ to the potential, see Eq. (fTSl i. 
in the s-wave phase shift fixing the scattering length ao to 
the same value. In view of the rather small g\ ~ 10"^ val- 
ues listed in table[T]we expect tiny corrections at even large 
energies kR(, ~ 10. This result not only shows a clear dom- 
inance of the leading dispersion coefficient Ce but also that 
scaling of VdW forces holds beyond naive dimensional es- 
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Fig. 5. Bound state spectrum for the VdW potential R\U(r) = 
-{R(,lrf[l + giRJrf] for 9 = (solid), g = 0.1 (long-dashed- 
dotted) ), 0.5 (dashed) and g = I (short-dashed-dotted). The hor- 
izontal line corresponds to a situation where a^lRf, = 0.336 and 
the intersections to y„ values with E„ = -yl/{2ji). 
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Fig. 6. s-wave phase shifts for the VdW potential R^U(r) = 
-{Re/rfll + giRilrf} for g = 0,0.1,0.5, 1. for a fixed sample 
value of the scattering length ao/R^ = 0.335. 

timates, kR(, ^ 1 . This raises the question on the usefulness 
of including higher order dispersion coefficients such as Cg 
and Cio since the region where they can distinctly be disen- 
tangled without entering the finite size regime is extremely 
narrow. This is also seen in the bound state spectrum dis- 
played in Fig. |5] since for yR(, < 10 a rather tiny energy 
shift is observed in the closest states to the continuum. 



5 Phenomenological potentials 

In our way of treating the renormalization of VdW forces, 
we need not specify the value of the scattering length, ao, 
since it exactly factors out in the expression for the phase 
shift (see Eq. dSTTl). However, to predict scattering phase 
shifts a particular value of ao must be used. It is interesting 
to analyze the results from a comparative perspective with 
the so called realistic inter-atomic potential models, which 
aim at a description through the entire range of distances. 
Thus, one can deduce from those the value of the scattering 
length. The advantage of such potentials is that they pro- 
vide a complete description of the interaction throughout 
the entire separation range. However, many of its features 
cannot be deduced accurately from first principles calcula- 
tions. In contrast, the long range part of the interaction has 
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a well accepted form in terms of a few parameters, say Ce, 
Cs, Cio which, in principle, can be evaluated from ab initio 
atomic structure electronic wave functions. 

Most modern inter-atomic potentials include the asymp- 
totic Van der Waals long distance behaviour. They are gen- 
erally written as the sum of a long range dispersive term 
and a short distance term with a core which reflects the 
impenetrability of two atoms. To keep the discussion as 
simple as possible in terms of the number of parameters, 
and for illustration purposes, we will analyze the venera- 
ble Lennard-Jones potential, which we write as 



2fiVuL = Uu(r) = 




where we have chosen to scale the potential in the long 
distance VdW units R = R(, = (IjiCe)* . The value of the 
dimensionless constant g determines the classical turning 
point UuigRe) - 0. In these form the minimum of the 
potential is at rmin -l^Rf, and - -l/(4gr 7?^. A rele- 
vant dimensionless parameter is the total number of bound 
states, which obviously increases as the repulsive term is 
shifted towards the origin. Within a WKB approximation 
the number of bound states is given by 

1 r°° 0.1339 

A^WKB = - dr^-Uuir) = ^ , (38) 

where a - gR(, is the zero energy classical turning point. 
For^ = 0.0365 we get A^WKB = 100, for g = 0.0517491 we 
get A^wKB - 50 and for g = 0.1 15715 we get A^wkb = 10. 

Using the VdW scaled units one can predict the scat- 
tering length ao from the LJ potential and the phase shifts 
unambiguously for any value of g. The result is displayed 
in Fig. |2] According to Levinson's theorem any time the 
scattering length jumps from -oo to -i-oo a new bound state 
dives from the continuum into the negative energy spec- 
trum. In Fig. |7] we display the effective range as given by 
Eq. ( [TT] ). The divergent values of ro correspond, accord- 
ing to our low energy theorem, to points where the scatter- 
ing length goes through zero. The very strong sensitivity 
to the precise location of the position of the core is man- 
ifest. In Fig. |2] we compare the universal and renormal- 
ized VdW eff'ective range formula with the actual values 
deduced from the Lennard-Jones potential for a different 
number of bound states = 1,5, 15 (plus or minus one). 
The minima in the curve corresponds to a value of the scat- 
tering length passing from -oo to +oo which corresponds 
to entering a new bound state in the spectrum. As one sees 
in the LJ case there is a multivalued function reflecting the 
multiple branches already observed in Fig. [T] The rather 
universal pattern of this figure is striking because it ex- 
plicitly shows that to very small uncertainties the value of 
the scattering length largely determines the value of the ef- 
fective range, regardless on the precise number of bound 
states. The only remarkable exception corresponds to the 
case with no bound states where the zero energy turning 
point is located at increasingly large distances. 

In a sense these correspond to almost low energy iden- 
tical situations, where additional bound states are hosted 
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Fig. 7. (Upper panel) The scattering length oq and effective range 
ro of the Lennard-Jones potential as a function of the zero en- 
ergy turning point gR(, (in VdW units). (Lower panel) The ef- 
fective range ro of the Lennard-Jones model as a function of the 
inverse scattering length I/q-q compared to the universal renor- 
malized VdW formula (in VdW units) for different number of 
bound states N = 1,5, 15. 



by the potential as the short range repulsion is displaced 
towards the origin. Moreover, one expects that the largest 
discrepancies from renormalized VdW and LJ should take 
place in the case of large scattering lengths, since a larger 
sensitivity to short distances is displayed in such a case. 
Actually, this is what happens. We note that there is a rela- 
tively constant shift between r^^ and r^ ''^ of about 0.5-1 .0 
for Iq-oI > IQR^. However, the renormalized theory works 
better the larger the number of bound states and also for 
small scattering lengths, since as we discussed in the pre- 
vious section for ao s> R(, large distances dominate. A 
relevant question is whether renormalization theory as ex- 
plained above can account for the behaviour of the phase 
shifts in a energy region where the De Broglie wavelength 
is larger than the short distance kR(,g <k 1 provided the 
scattering length ao is the same. The result is shown in 
Fig.[8]for a variety of values of g which cover several cases 
with large and small scattering lengths as compared to the 
VdW radius. We see the anticipated similarity despite the 
fact that both potentials are completely different at short 
distances. Actually, the phase-shifts are indistinguishable 
for kR(y <K 1, but they go hand in hand far beyond this 
expected region; what matters is kR(,g <s; 1 . On the other 
hand, the ERE given by Eq. ( fTSI ) truncated to second order, 
i.e. taking U2 = only works for kR(, <s: 1. In a sense this 
is equivalent to "seeing" the Van der Waals force in a scat- 
tering experiment. The point of renormaUzation theory is 
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Fig. 8. The Lennard- Jones s-wave phase shifts So(k) (in radians) for different number of bound states N = 0, 1,3, 5 as a function of 
the wave number (points) compared with the corresponding VdW renormalized ones having the same scattering length ao (solid). We 
compare also with the effective range expansion truncated to second order k cot do(k) = - l/ag + rgk^/l (dots). 



that it yields model independent results and hence any dis- 
crepancy with data can be clearly attributed to a deficient 
incorporation of the long distance physics. 



6 The Effective Field Theory and its limits 

At very low energies the interaction between neutral atoms 
can be handled by an eff'ective range expansion, Eq. (fTSl l. 
where the long range character of the VdW potential be- 
comes manifest in the third term of the expansion. How- 
ever, if only the first two terms are retained 



1 1 , 

kcot6o(k) 1 — ri)k . 

ao 2 



(39) 



there arises the interesting possibility of universally repre- 
senting long range forces on equal footing with short range 
interactions. We may judge the quality of such an appeal- 
ing approximation by comparing Fig. |8]the result of the 
renormalized VdW theory with the ERE. As we see the 
expansion truncated to second order breaks down at low 
energies, namely kR(, <K 1, as expected. 

Under these very restrictive conditions the problem can 
be advantageously treated by using EFT methods, which 



are based on the compelling idea that at such long wave- 
lengths atoms behave as elementary structureless particles. 
This point of view has been stressed recently (see e.g. [ 14 15]) 
with particular fruitful predictions in the three-body prob- 
lem 1 50 16 1 where Efimov states have been predicted. This 
is usually done by considering the Lagrangian density 



id, + 



2m 



iA(x) 



tfr\x)4^\x')V(x - x')4r(x')tfr(x) , (40) 



where i^(x) = ifr(x, t) are space-time canonically quantized 
fields with Fermi or Boson statistics depending upon the 
spin nature of the atom as a whole. The irreducible two- 
point function is the potential V{x) which is taken as 

<k'|V|k) = J d^xe'^^-^'^''V(x)^Co + C2(k'^ + k^)(4-l) 

Usually the coefficients Co and C2 are assumed to be com- 
pletely arbitrary. The corresponding Lippmann-Schwinger 
equation can be solved requiring introducing a momentum 
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space cut-off A (see e.g. Ref. lISTl and references therein) 



J_ _ 10(C2MA^ - 3f 
~ 9Mn(-C^MA^ + 5Co) 



2A 
n 



(42) 



1 50C2f3 -HCaMA^) (6 + CaMA^) ? 

^ ' ^ - + -r .(43) 

TlA 



TlTi[-iCQ + C\A^M) 



This leads for any cut-off A to the mapping (oro, ro) — » 
(Co, C2). Eliminating Co and C2 in favour of co and ro, 
the phase shift becomes 



2A 



(n - lAa^f 



TTffo 1A{ti - 2Aao) + Q'op^(ro;rA - 4) 



2A p A + p 

+ - log 

7T n A - p 



(44) 



which is a real quantity for p < A, meaning that two-body 
unitarity is fulfilled. However, one has complex solutions 
for Co and C2 if 

alronA^ - l6alA^ + UaonA - Stt^ < . (45) 

which means that the effective Lagrangian, Eq. ( l40l l. be- 
comes non hermitian -C^ix) + -C{x). On the other hand, 
it is well known that three-body unitarity rests on off-shell 
two-body unitarity ll52l . a condition which cannot be met if 
the interaction is not hermitian, mainly because Schwartz's 
principle. So the lesson is that while a complex two-body 
potential may fulfill on-shell two-body unitarity, a viola- 
tion of three body unitarity is still possible 1.5 1 j . 

Generally, Eq. (|45] | imposes a limit on the maximum 
value of the cut-off A, but to find it we need to know both 
cKo and ro. For our case of VdW interactions we have seen 
that the universal formula for ro in terms of ao, Eq. ( |32] | 
works extremely well if we judge by Fig. |2] Thus, if we 
merge Eq. (l45l l equal to null and Eq. ( |32] | we obtain a 
boundary in the (A, ao) plane which is suitably represented 
in Fig.|9] Essentially and up to minor variations the mean- 
ing is that the cut-off A cannot exceed the VdW wave num- 
ber, A < 7t/{2R(,). This issue is relevant for three-body cal- 
culations |53 54| as addressed recently [55|. 

The analysis of the problem in coordinate space assum- 
ing an effective local and energy independent long distance 
dynamics and an energy dependent boundary condition at 
short distances have been analyzed in Ref. [26 1 for s-waves 
and in Ref. [561 in the three-dimensional case. It is found 
that the locality condition for an s-wave implies 



d 



m;,s('-c) 



Uk.&{rc) 



< 



(46) 



where Uk s(r) is the wave function for r < r^. Note the 
resemblance with Eq. (fTOl i. If there would be no interaction 
for r > Tc then we have 

Uk{r) = sin(^r + 6i)(k)) , r > rc , (47) 

which can be matched to the inner r < r^. region yielding 

d 



dk^ 



[kcot{kr,+6o(k))] < 0, 



(48) 
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Fig. 9. Boundaries on the cut-off for Van der Waals systems with 
reduced potential IjiVir) = -R'*/r^ for which causahty in coordi- 
nate space (left panel) and off-shell unitarity in momentum space 
(right panel) is preserved. 



this is equivalent to Wigner's causality condition fSTl, as 
noted in |58|, and combined to the ERE, Eq. ( [39] l, provides 
a constraint on the effective range 



ro/(2r,) < 1 - rjao + r^J(3al) . 



(49) 



If there is no potential for r > rc and take ro = ryaw 
we get a universal lower limit for r^. . A similar conclusion 
has been presented recently [59|. The conditions featuring 
causality in coordinate space as well as off-shell unitarity 
in momentum space are depicted in Fig. |9] and suggest 
that modelling a finite range of VdW forces by an effective 
Lagrangian such as Eq. ( l40l l requires assuming a cut-off 
distance larger than the VdW length. 



7 Conclusions 

Renormalization ideas can profitably be exploited in con- 
junction with the superposition principle of boundary con- 
ditions in the description of model independent and univer- 
sal features of the VdW force. Our main points are 

- Van der Waals interactions between neutral atoms obey 
scaling rules which allow to determine the scattering 
and binding properties universally. They are well satis- 
fied phenomenologically and extend much beyond low 
energy approximations such as the effective range ex- 
pansion. 

- There is a clear dominance of the leading C(, contri- 
bution in a rather wide energy range. The range where 
higher order corrections due to Cg or Cio provide a dis- 
tinct correction yet the finite size effects can still be 
neglected is extremely narrow or inexistent. 

- Van der Waals potentials can be represented by short 
distance contact interactions under restrictive condi- 
tions based on causality and off-shell unitarity which 
are independent on the value of the scattering length. It 
is inconsistent to model VdW forces assuming a cut-off 
distance smaller than the VdW length. 

We thank M. Pavon Valderrama and R. Gonzalez Ferez 
for discussions. Work supported by Spanish DGI and FEDER 



EPJ Web of Conferences 



funds with grant FIS2008-01 143, Junta de Andalucia grant 
FQM-225-05, EU Integrated Infrastructure Initiative Hadron 
Physics Project conti'act RII3-CT-2004-506078. 



References 



1. J.S. Rowlinson, Cohesion (Cambridge University 
Press, 2005). 

2. E London, Zeitschrift fur Physik 63 (1930) 245. 

3. H.B.G. Casimir and D. Polder, Phys. Rev. 73 (1948) 
360. 

4. G. Feinberg and J. Sucher, Phys. Rev. A2 (1970) 2395. 

5. G. Eeinberg, J. Sucher and C.K. Au, Phys. Rept. 180 
(1989)83. 

6. K.M. Case, Phys. Rev. 80 (1950) 797. 

7. W. Erank, D.J. Land and R.M. Spector, Rev. Mod. 
Phys. 43 (1971)36. 

8. J. Weiner, Cold and Ultracold Collisions in Quantum 
Microscopic and Mesoscopic Systems (Cambridge 
University Press, 2004). 

9. G.E Gribakin and V.V. Elambaum, Phys. Rev. A 48 
(1993)546. 

10. B. Gao, Phys. Rev. A58 (1998) 1728. 

11. V.V. Elambaum, G.E. Gribakin and C. Harabati, Phys. 
Rev. A 59 (1999) 1998. 

12. B. Gao, Phys. Rev. A 80 (2009) 012702. 

13. A. Calle Cordon and E. Ruiz Arriola, (2009), 
09 1 2 . 1 7 1 4 [cond-mat .other] . 

14. E. Braaten and H.W. Hammer, Phys. Rept. 428 (2006) 
259, cond-mat/0410417 

15. E. Braaten, M. Kusunoki and D. Zhang, Annals Phys. 
323 (2008) 1770, 0709.0499. 

16. L. Platter, Eew Body Syst. 46 (2009) 139, 0904.2227. 

17. I.N. Levine, Quantum Chemistry (Allyn and Bacon, 
1974). 

18. W. Heitler and E. London, Zeitschrift fur Physik 44 
(1927)455. 

19. W. Kolos and C.C. Roothaan, Reviews of Modern 
Physics 32(1960)219. 

20. R. Eisenschitz and E. London, Zeitschrift fur Physik 
60(1930) 491. 

21. A. Calle Cordon, Renormalizacion de Interacciones 
Atomicas mediante condiciones de contorno (Master 
Thesis (University of Granada) December, 2007). 

22. J. Mitroy and M.W.J. Bromley, Phys. Rev. A 71 (2005) 
032709, arXiv:physics/041 1 172 

23. E.H. Lieb and W.E. Thirring, Phys. Rev. A 34 (1986) 
40. 

24. K. Cahill and VA. Parsegian, eprint 
ai-Xiv:q-bio/03 12005 (2003), arXiv:q-bio/03 12005 

25. A.M. Lane and R.G. Thomas, Rev. Mod. Phys. 30 
(1958)257. 

26. M. Pavon Valderrama and E.R. Arriola, Phys. Rev. 
C74 (2006) 054001, nucl-th/0506047 

27. M. Pavon Valderrama and E.R. Arriola, Annals Phys. 
323 (2008) 1037, 0705.2952. 

28. E. Ruiz Arriola, A. Calle Cordon and M. 
Pavon Valderrama, (2007), 0710.2770. 



29. S.V Khristenko, A.I. Maslov and VP Sheveiko, 
Molecules and Their Spectroscopic Properties 
(Springer Series on Atoms and Plasmas. Vol. 21. 
Springer- Verlag, 1998). 

30. Z. Yan et al., Phys. Rev. A 54 (1996) 2824. 

31. A. Derevianko, J.E Babb and A. Dalgarno, Phys. Rev. 
A 63 (2001) 052704, arXiv:physics/0102030 

32. S.G. Porsev and A. Derevianko, Jour. Chem. Phys. 
1 19 (2003) 844, |arXiv:physics/0303048 l 

33. M. Marinescu, D. Vrinceanu and H.R. Sadeghpour, 
Phys. Rev. A 58 (1998) 4259. 

34. S.G. Porsev and A. Derevianko, Soviet Journal of Ex- 
perimental and Theoretical Physics 102 (2006) 195. 

35. J. Werner et al.. Physical Review Letters 94 (2005) 
183201, arXiv:cond-mat/0412049 

36. B.R. Levy and J.B. Keller, Journal of Mathematical 
Physics 4 (1963) 54. 

37. R. Cote, E.J. Heller and A. Dalgarno, Phys. Rev. A 53 
(1996) 234. 

38. R. Cote and A. Dalgarno, Phys. Rev. A 50 (1994) 
4827. 

39. V.V. Elambaum, G.E Gribakin and C. Harabati, Phys. 
Rev. A 59 (1999) 1998. 

40. M. Marinescu, Phys. Rev. A 50 (1994) 3177. 

41. H. Ouerdane and M.J. Jamieson, European Physical 
Journal D 53 (2009) 27, 0802.1222. 

42. M.J. Jamieson et al.. Journal of Physics B Atomic 
Molecular Physics 40 (2007) 3497. 

43. M.J. Jamieson et al.. Journal of Physics B Atomic 
Molecular Physics 36 (2003) 1085. 

44. Z. Pavlovic et al., Phys. Rev. A 69 (2004) 030701, 
arXiv:physics/0309076 

45. M. Kemal Oztiirk and S. Ozgelik, ArXiv Physics e- 
prints (2004), arXiv:physics/0406027 

46. N. Koyama and J.C. Baird, Journal of the Physical 
Society of Japan 55 (1986) 801. 

47. M.J. Jamieson, A. Dalgarno and J.N. Yukich, Phys. 
Rev. A 46 (1992) 6956. 

48. A. Sen, S. Chakraborty and A.S. Ghosh, Europhysics 
Letters 76 (2006) 582. 

49. M.J. Jamieson, A. Dalgarno and M. Kimura, Phys. 
Rev. A 51 (1995) 2626. 

50. E. Braaten and H.W. Hammer, Annals Phys. 322 
(2007) 120, cond-mat/0612123 

51. D.R. Entem et al., Phys. Rev. C77 (2008) 044006, 
0709.2770. 

52. C. Lovelace, Phys. Rev. 135 (1964) B1225. 

53. L. Platter and D.R. Phillips, Eew Body Syst. 40 (2006) 
35, cond-mat/0604255 

54. L. Platter, C. Ji and D.R. Phillips, Phys. Rev. A79 
(2009) 022702, 0808.1230. 

55. D.L. Canham and H.W. Hammer, (2009), 0911.3238. 

56. A. Calle Cordon and E. Ruiz Amola, (2009), 
0905.4933[nucl-th]. 

57. E.P Wigner, Phys. Rev. 98 (1955) 145. 

58. D.R. Phillips and T.D. Cohen, Phys. Lett. B390 (1997) 
7,|nucl-th/9607048 

59. H.W. Hammer and D. Lee, Phys. Lett. B681 (2009) 
500, 0907.1763. 



